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ON Abstract 

00 . 

f^**) Here we show that a particular one-parameter generalization of the exponential function is suitable to unify most of the popular 
CO one-species discrete population dynamics models into a simple formula. A physical interpretation is given to this new introduced 
parameter in the context of the continuous Richards model, which remains valid for the discrete case. From the discretization of 
f***) the continuous Richards' model (generalization of the Gompertz and Verhuslt models), one obtains a generalized logistic map and 
OO we briefly study its properties. Notice, however that the physical interpretation for the introduced parameter persists valid for the 
discrete case. Next, we generalize the (scramble competition) #-Ricker discrete model and analytically calculate the fixed points 
as well as their stability. In contrast to previous generalizations, from the generalized 0-Ricker model one is able to retrieve either 
. _ scramble or contest models. 

^ Key words: Complex Systems, Population dynamics (ecology), Nonlinear dynamics 
- PACS: 89.75.-k, 87.23.-n, 87.23.Cc, 05.45.-a 



1. Introduction 

Recently, the generalizations of the logarithmic and ex- 
ponential functions have attracted the attention of resear- 
chers. One-parameter logarithmic and exponential func- 
tions have been proposed in the context of non-extensive 
statistical mechanics |H2lol 115] . relativistic statistical me- 
chanics |6l7j and quantum group theory [8] . Two and three- 
parameter generalization of these functions have also been 
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proposed [9110111] . These generalizations are in current use 
in a wide range of disciplines since they permit the gener- 
alization of special functions: hyperbolic and trigonomet- 
ric pj], Gaussian/Cauchy probability distribution func- 
tion [T5] etc. Also, they permit the description of several 
complex systems |14I15I16I17I18 19 , for instance in gener- 
alizing the stretched exponential function . 

As mentioned above, the one-parameter generalizations 
of the logarithm and exponential functions are not univo- 
quous. The q- logarithm function hiq(x) is defined as the 
value of the area underneath the non-symmetric hyperbole, 
f~(t) = l/t 1 -*, in the interval t G [1, x] [21]: 



f dt x q — 1 

J t 1 1 q'^q q 
1 

This function is noi the ordinary logarithmic function in 
the basis namely [log^(a^)], but a generalization of the 
natural logarithmic function definition, which is recovered 
for q — 0. The area is negative for < x < 1, it vanishes 
for x = 1 and it is positive for x > 1, independently of the 
g values. 

Given the area a; underneath the curve fq(t), fort 6 [0,y], 
the upper limit ?/ is the generalized q-exponential function: 
y = e q (x). This is the inverse function of the q- logarithmic 
eg[hig(a;)] = x = hiq[eq(x)] and it is given by: 

!0 for qx < — 1 

, m- ■ ( 2 ) 

lim(l + gV) 1/9 forgx>-l 
q'^q 

This is a non-negative function e,j(x) > 0, with eq(0) = 1, 
for any <j. For q — > ±oo, one has that e-oofx) = 1, for x < 
and 600(2:) = 1, for x > 0. Notice that letting a; = 1 one 
has generalized the Euler's number: 

e q (l) = (l + q) 1 /*. (3) 

Instead of using the standard entropic index q in Eqs. (flj) 
and (jnj, we have adopted the notation g = 1 — q. The latter 
notation permits us to write simple relations as: hiq(x) = 
— hx—q{x) or e_q(— x) — l/eq(x), bringing the inversion 
point around q — 0. These relations lead to simpler expres- 
sions in population dynamics problems [25) and the general- 
ized stretched exponential function [20] contexts. Also, they 
simplify the generalized sum and product operators |21j . 
where a link to the aritmethical and geometrical averages 
of the generalized functions is established. 

This logarithm generalization, as shown in Ref. (22] p. 
83] , is the one of non-extensive statistical mechanics [2] . It 
turns out to be precisely the form proposed by Montroll 
and Badger [23] to unify the Verhulst (q = 1) and Gom- 
pertz (q = 0) one-species population dynamics model. 
The g-logarithm leads exactly to the Richards' growth 
model [23125] : 

dlnpjt) 

— - — =-Khigp{t), (4) 

where p(t) — N(i)/Noo, N(i) is the population size at 
time t, Nov is the carrying capacity and k is the in- 
trinsic growth rate. The solution of Eq. (j4]) is the q- 
generalized logistic equation p(t) — l/e q -{\nq(pQ 1 )e~ Kt ] — 
e_ 5 [-ln ? (po 1 )e-' ,t ] = e^ q [ln^ q (p ) e - Kt ]. 

The competition among cells drive to replicate and 
inhibitory interactions, that are modeled by long range 
interaction among these cells. These interactions furnish 
an interesting microscopic mechanism to obtain Richards' 
model [26127] . The long range interaction is dependent on 
the distance r between two cells as a power law r 7 . These 
cells have a fractal structure characterized by a fractal 
dimension Df. 



Here we call the attention to Eq. (7) of Ref. [25], 
namely h(t) = n(t){(G) - JI[n(t)}}, where I(n(t)) = 
w{[ J D / n(t)/u;] 1 -T'/ D / - 1} /[D f (l - 7/-D/)]. Here, w is a 
constant related to geometry of the problem, (G) is the 
mean intrinsic replication rate of the cells and J is the 
interaction factor. Using Eq. ([1]), one can rewrite it sim- 
ply as: dlnn(*)/dt = (G)/n(t) - Jw\nq[Dfn(t)/u]/D f . 
Calling, p — Dfri/u), n — Juj/Df and q = 1 — j/Df, this 
equation is the Richard's model [Eq. with an effort rate 
(G)/n(t). In this context the parameter q acquires a phys- 
ical meaning related to the interaction range 7 and fractal 
dimension of the cellular structure Df. If the interaction 
does not depend on the distance, 7 = 0, and it implies 
that q = 1 . This physical interpretation of q has only been 
possible due to Richards' model underlying microscopic 
description. 

Introduced by Nicholson in 1954 28J, scramble and con- 
test are types of intraspecific competition models that differ 
between themselves in the way that limited resources are 
shared among individuals. In scramble competition, the re- 
source is equally shared among the individuals of the pop- 
ulation as long as it is available. In this case, there is a 
critical population size N c , above which, the amount of re- 
source is not enough to assure population survival. In the 
contest competition, stronger individuals get the amount of 
resources they need to survive. If there is enough resources 
to all individuals, population grows, otherwise, only the 
strongest individuals survive (strong hierarchy), and the 
population maintains itself stable with size . 

From experimental data, it is known that other than the 
important parameter k (and sometimes -ZVqo), additional 
parameters in more complex models are needed to adjust 
the model to the given population. One of the most general 
discrete model is the #-Ricker model [29J30J. This model 
describes well scramble competition models but it is unable 
to put into a unique formulation the contest competition 
models such as Hassel model [28] , Beverton-Holt model [31] 
and Maynard-Smith-Slatkin model [32] ■ 

Our main purpose is to show that Eq. ^ is suitable to 
unify most of the known discrete growth models into a sim- 
ple formula. This is done in the following way. In Sec. [21 we 
show that the Richards' model [Eq. {3}], which has an un- 
derlying microscopic model, has a physical interpretation 
to the parameter q, and its discretization leads to a gener- 
alized logistic map. We briefly study the properties of this 
map and show that some features of it (fixed points, cycles 
etc.) are given in terms of the g-exponential function. Cu- 
riously, the map attractor can be suitably written in terms 
of q-exponentials, even in the logistic case. In Sec. [3] us- 
ing the ^-exponential function, we generalize the 0-Ricker 
model and analytically calculate the model fixed points, 
as well as their stability. In Sec. (U we consider the gen- 
eralized Skellam model. These generalizations allow us to 
recover most of the well-known scramble/contest competi- 
tion models. Final remarks are presented in Sec[SJ 



2. Discretization of the Richards' model 



To discretize Eq. (gj), call (p 



Pi) /At 



1)/?, 
to: 



p~ = l + kAt/q&ndxi = p i [(p q - 



1)/Pd 9 



x i+l = p'q X iQ - X I) = ~Pq x i foq(Xi) , 



which leads 



(5) 



where pq = qp'^. We notice that q keeps its physical inter- 
pretation of the continuous model. 

In Eq. if q = 1 and pi = p[ — 4a, with a S [0, 1], one 
obtains the logistic map, iCj+i = 4aa;,(l — CEj), which is the 
classical example of a dynamic system obtained from the 
discretization of the Verhulst model. Although simple, this 
map presents a extremely rich behavior, universal period 
duplication, chaos etc. [33] . 

Let us digress considering the Feigenbaum's map [34] : 
j/j+i = 1 — py1 +1 , with q > 0, < p < 2 and — 1 < 
Hi < 1- Firstly, let us consider the particular case q = 1. 
If one writes yi = yi — b, with being a constant, then: 
= — pb 2 + b + 1 + 2/i6j/i[i — yi/(2b)]. Imposing that 
-pb 2 + b 
calling Xi - 
4a = 



1 = leads to b± = (1 T a/ 1 + 4 /-0/( 2 A*) and 
yi/(2b), one obtains the logistic map with pi = 



= 4a = 1 + t/1 + 4p, so that < pi < 4. One can eas- 
ily relate the control parameter of these two maps, mak- 
ing the maps equivalent. For arbitrary values of q, there is 
not a general closed analytical form to expand — b\^ +1 
and one cannot simply transform the control parameters of 
Eq. to the Feigenbaum's map. Here, in general, these 
two maps are not equivalent. It would be then interesting, 
to study the sensitivity of Eq. §5§ with respect to initial 
conditions as it has been extensively studied in the Feigen- 
baum's map [35136137138] . 

Returning to Eq. ([5]), in the domain < x < 1, f(xi) = 
— pqXi hiq(xi) > (non-negative), for pq > 0. Since eq(x) 
is real only for qx > — 1, / is real only for q > — 1. The 
maximum value of the function is 

Pq 



f = m 



which occurs at 
1 



x 



e*(l) 



(6) 



(7) 



i. e., the inverse of the generalized Euler's number eq(l) 
[Eq. (|3J)]. For the generalized logistic map, < x < 1, so 
that < / < 1, it leads to the following domain for the 
control parameter < pq < p m ax'- 



Pmax = eg(l)ei(g) = (1 
The map fixed points [x* 

xl=eq{-llPq) . 



f(x*)} are 



(8) 

(9) 
(10) 



The fixed point x\ is stable for < pq < q and x\ is stable 
for q < p q < p P d: where 



Notice the presence of the g-exponentials in the description 
of the attractors, even for the logistic map 5=1. 

The generalized logistic map also presents the rich be- 
havior of the logistic map as depicted by the bifurcation di- 
agram of Fig.[T] The inset of Fig. [1] displays the Lyapunov 
exponents as function of the central parameter p^. 
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Fig. 1. Bifurcation diagram of Eq. (0 for q = 2, where we see that the 
period doubling start at p = (q + 2)/q = 2 [Eq. 111H and the chaotic 
phase finishes at p ma x = e 2 (l)e 1 (2)/q = V27/2 a 2.6 [Eq. ©]. 
Inset: Lyapunov exponents as function of pq for q = 2. 

In Fig.[5]we have scaled the axis to pq[—(q+2) /q]/ (p m axq), 
where p ma x is given by Eq. ([5]) and we plotted the bifur- 
cation diagram for q — 1/10, 1 and 10. We see that the 
diagrams display the same structure but each one has 
its own scaling parameters. The role of increasing q is 
to lift the bifurcation diagram to relatively anticipating 
the chaotic phase. The period doubling region start at 
x*(q) = e q [-l/(q + 2) = [1 - 1/(1 + 2/q)} 1 /^ so that for 
x{{1/10) = (20/21) 10 w 0.61, x* 2 (l) = 2/3 w 0.67 and 
x*(10) = (1/6) 1 / 10 m 0.84. 
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Fig. 2. Bifurcation diagram of Eq. {5} for q = 1/2, q = 1 (logistic) 
and q = 2. The fixed points are given by Eqs. (J9j and 1101 . 

When = eq(l)ei(q), then a;^ e (0,1). In Fig. [3] we 
show the histograms of the distribution of the variable Xi. 



We see that as q increases, the histograms have the same 
shape as the logistic histogram has, but it is crooked in the 
counter clock sense around x = 1/2. 




Fig. 3. Histograms of x for p<j = eq{l)e\{q) in Eq. J5} and q = 1/10, 
q = 1 (logistic) and q = 10. There is not a drastic change with 
respect to the logistic map. Only on the corners one is able to see 
the difference as shown by the insets. 



3. The generalized 0-Ricker model 

The 9 -Richer model [29.30J is given by: 



T . M _ T . p r[l-(xi/K) e ] 



(12) 



where 6 > 0. 

Notice that x — r 1 ! 8 xj k is the relevant variable, where 
K\ — e r > 0. In this way Eq. (|12p can be simply writ- 
ten as ii+i = k\Xie~ Xi . For 9 = 1, one finds the standard 
Richer model [40] . For arbitrary 9, expanding the exponen- 
tial to the first order one obtains the generalized logistic 
map [Eq. ([SJ] which becomes the logistic map, for 9 = 1. 
The 0-Ricker, Ricker and quadratic models are all scramble 
competion models. 

If one switches the exponential function for the ^-gene- 
ralized exponential in Eq. p2[) . one gets the generalized 9- 
Richer model: 



V K 



K\Xi 



l + qr^Y 



1/9 



(13) 



To obtain standard notation, write c = 1/q and &2 = 
r/(kc), so that Xi+\ = k\Xi/{l + k2Xi) c [41] . 

The generalized model with 9 = 1, leads to the Hassel 
model [55] , which can be a scramble or contest competition 
model. One well-known contest competition model is the 
Beverton-Holt model [21] , which is obtained taking q = c = 
1. For q = 0, one recovers the Ricker model and for q = — 1, 
one recovers the logistic model. 

It is interesting to mention that the Beverton-Holt 
model [31] is one of the few models that have the time evolu- 
tion explicitly written: Xi = k\xq/\1 + (1 — k\)xq/(1 — Ki)]. 



From this equation, one sees that £i(>i) = ,for Ki < 1 
and »i(>i) = Kx — 1 for k\ > 1. 

Using arbitrary values of 9 in Eq. (fT3|) . for q = one 
recovers the 9- Ricker model, and for q = 1, the Maynard- 
Smith-Slathin model |32] is recovered. The latter is a scram- 
ble/contest competition model. For q = — 1, one recovers 
the generalized logistic map. The trivial linear model is re- 
trieved for q —> — oo. 

In terms of the relevant variable x, Eq. (| 13[) is rewritten 
as: 



Xi+l 



i ~0\ 

K\XiC — q\ j , 



(14) 



where Xi > and we stress that the important parame- 
ters are K\ > 0, q and 6 > 0. Eq. (fF4"]) is suitable for data 
analysis and the most usual known discrete growth models 
are recovered with the judicious choice of the q and 9 pa- 
rameters as it shown in Table [1] Some typical bifurcation 
diagrams of Eq. (Ti~4"]) are displayed in Fig. [3J 
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<J 
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competition 
type 


Linear 


— oo 


> 




Logistic 


-1 


i 
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Generalized Logistic 


-1 
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s 


Ricker 





1 


s 


6>-R,icker 





> 


s 


Hassel 


* 


1 


s (q < 1/2) or c (q > 1/2) 


Maynard-Smith-Slatkin 


1 


> 


s (0 > 2 ) or c (6 < 2) 


Beverton-Holt 


1 


1 


c 



Table 1 

Summary of the parameters to obtain discrete growth models from 
Eq. (Tl4]l. In the competition type column, s and c stand for scramble 
and contest models, respectively. The symbol * stands for arbitrary 
values. 

Now, let us obtain some analytical results for the map of 
Eq. (|14p . which we write as ii+i = f g t r {xi), with 



fgtr(x) = KlXe^q(~X ) = 



K\X 



(1 + qx 6 ) 1 ^ 



(15) 



The i domain is unbounded (x > 0), for q > 0. However, 
for q < 0, fgtr(x) = kix(1 — Iqlx 6 ) 1 /^ and the i-domain 
is bounded to the interval: < x < x m , with 

' (16) 



*" m ~ (-q)V e ' 

so that for |g| < 1, x m > 1; for |g| = 1, with x m = 1 
(for q = —1, it is the generalized logistic case [Eq. $5$i in 
9 instead of q] and for q = 1, the Maynard-Smith-Slatkin 
model) and for |<7| > 1, x m < 1. 

The derivative of f g t r with respect to x is: fg tr (x) = 
K 1 [l + (q-9)x e ]/(l + qx e ) 1 + 1 ^. Imposing f gtr (x max ) = 0, 
one obtains the maximum of Eq. (|15|) . 



fgtr — fgtr (Xmax ) 



nix m e q (-l/9) 
ee(-l/q) 



(17) 




(a) 



ee(-l/q) 



e e (-l/g)e_ 9 (l/0) 



(19) 



From Eq. (US]), one sees that for 9 < q, f g t r (x) does 
not have a hump, it is simply a monotonically increasing 
lunction of x, which characterizes the contest models. At 
the critical 9 = q value, the function f g tr(x) starts to have a 
maximum value at infinity. For 9 > q, the function f g tr(x) 
has a hump, with maximum value at x ma x [Eq. (|18p ] such 
that as 9 — ► oo then x max — ► 1. 

The map fixed points [x* = f g tr(x*)} are: 





(c) 



300 



Fig. 4. Typical bifurcation diagrams of Eq. 1140 . (a) q = 1/2 and 
= 2, (a) q = 1 and = 3 and (a) q = 2 and = 5. 

at 



'-ff)V« ee(-l/q) 



(18) 



The control parameter Hi is unbounded («i > 0), for q > 0, 
but for q < 0, since / g t,- < x rn , it belongs to the interval 
< K% < K m , where: 



x{ = Q 

x* 2 = K-(ki)] 1 / 9 > 



(20) 
(21) 



These fixed points are show as function of K\ in Fig. 




Fig. 5. Fixed points, given by Eqs. I|20| l and l|21| l. as function of fei. 
Instability regions are represented by dashed lines. The value of re m 
is given by Eq. H19I I. 

The fixed point x\ represents the species extinction and 
is stable, for < K\ < 1. Both fixed points and are 
marginal, for Ki = 1. For 1 < Ki < e- g (2/9), x\ becomes 
unstable and x\ > is stable and represents the species 
survival. For m = e- g (2/9), x* 2 becomes unstable and as 
K\ increases, a stable cycle-2 appears. For q < 0, as k\ in- 
creases further, the cycle-2 becomes unstable at some value 
of K\ giving rise to a route to chaos as in the logistic map, 
via period doubling. Nevertheless, for q > 0, several sce- 
narios may take place. Even though q > and 9 > q, if q < 
9 < 2q, the map f g t r has a hump, but it is not thin enough 
to produce periods greater than unity. In this case, f g t r 
produces only the two fixed points i* and xi,-, which char- 
acterize the context models. Nevertheless, scramble mod- 
els (9 > 2q) have maps with a hump thin enough to pro- 
duce stable cycles with period greater than unity. Thus, in 
scramble models one has more complex scenarios such as 
period doubling, as a route to chaos, as K\ > e- g {2/9). We 
have not being able to obtain analytically the behavior of 
the system q > and 9 > 2q. The 9 xq diagram is depicted 
in Fig. [6] 

For q = 6 = 1, one retrieves the Beverton-Holt model, 
with the fixed points x\ — and x\ — hii(Ki) = K\ — 1. 
For q, one retrieves the generalized logistic map. 



r. 




-1 1 



x m = {-^\(-V(fj 1 - *,= » 

Fig. 6. Diagram distinguishing the several types of behavior of the 
generalized 9-Ricker model [Eq. I|14| l] . For q < 0, one has the map 
= KiXi(l — I^IS?) 1 ' and < x < x m , with x m given by 
Eq. I|16j l, This system represents scramble models and achieves chaos, 
period doubling. For q > 0, x > and the behavior of the system is 
driving by the 9 values. For 9 > 2q, we do no have analytical results. 
Numerical simulation indicate regions with finite order of period 
doubling. For 9 < 2q, the system represent context models, with 
only two fixed points given by Eqs. (1206 and 121> . The map presents 
a hump, for 9 > 2q, but for 9 < q, the map is a monotonically 
increasing function. 

4. Generalized Skellam model 



structure. From the discretization of the Richard's model, 
we have obtained a generalization for the logistic map and 
briefly studied its properties. An interesting generalization 
is the one of 6*-Ricker model, which allows to have several 
scramble or contest competition discrete growth models as 
particular cases. Equation (fT4|) allows the use of softwares 
to fit data to find the most suitable known model through- 
out the optimum choice of q and 9. Furthermore, one can 
also generalize the Skellam contest model. Only a few spe- 
cific models mentioned in Ref. [41] are not retrieved from 
our generalization. Actually, we propose a general proce- 
dure where we do not necessarily need to be tied to a spe- 
cific model, since one can have arbitrary values of q and 9. 



Acknowledgments The authors thank C. A. S. Tergariol 
for fruitful discussions. ASM acknowledges the Brazil- 
ian agency CNPq (303990/2007-4 and 476862/2007-8) or 
support. RSG also acknowledges CNPq (140420/2007-0) 
for support. ALE acknowledges CNPq for the fellowship 
and FAPESP and MCT/CNPq Fundo Setorial de Infra- 
Estrutura (2006/60333-0). 



All the contest competition models generalized by 
Eq. (|13[) are power-law-like models for g / 0. However, 
the Skellam contest model cannot be obtained from this 
approach. It is the complement of an exponential decay 
= k(1 — e~ rXi ) [42]. Nevertheless, it is interesting to 
replace the exponential function to the ^-exponential in 
this model: Xi+i — k[l — e_,j(— rxi)) and write x = rx and 
k = rk, which leads to: 

Xi+i = k [1 - e-g(-Xi)] . (22) 

For q — oo, Eq. (|22|) leads to the constant model, for 
q = — 1, the trivial linear growth is found. If q = 0, one 
recovers the Skellam model and finally, q = 1 leads to the 
Beverton-Holt contest model (see Table [2]). 



Model 


q 


constant 


—00 


linear 


-1 


Skellam 





Beverton-Holt 


1 



Table 2 

Summary of the parameters to obtain contest competition discrete 
growth models from Eq. (I22t . 

5. Conclusion 

We have shown that the ^-generalization of the exponen- 
tial function is suitable to describe discrete growth mod- 
els. The q parameter is related to the range of a repulsive 
potential and the dimensionality of the fractal underlying 
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